Morphological instability of core-shell metallic nanoparticles 
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^ ; Abstract 

o ' Bimetallic nanopaiticles (often known as nanoalloys) with core-shell aiTangement are of special interest 



in several applications, such as in optics, catalysis, magnetism and biomedicine. Despite wide interest in 
applications, the physical factors stabilizing the structures of these nanoparticles are still unclear to a great 
extent, especially for what concerns the relationship between geometric structure and chemical ordering 
(-| ■ pattern. Here global-optimization searches are performed in order to single out the most stable chemical 



ordering patterns corresponding to the most important geometric structures, for a series of weakly miscible 
systems, including AgCu, AgNi, AgCo and AuCo. The calculations show that (i) the overall geometric 
structure of the nanoalloy and the shape and placement of its inner core are strictly correlated; (ii) centered 
cores can be obtained in icosahedral nanoparticles but not in crystalline or decahedral ones, in which asym- 
metric quasi-Janus morphologies form; (iii) in icosahedral nanoparticles, when the core exceeds a critical 
size, a new type of morphological instability develops, making the core asymmetric and extending it to- 
wards the nanoparticle surface; (iv) multi-center patterns can be obtained in polyicosahedral nanoalloys. 
Analogies and differences between the instability of the core in icosahedral nanoalloys and the Stranski- 



> 

■^ \ Krastanov instability occurring in thin-film growth are discussed. All these issues are crucial for designing 



strategies to achieve effective coatings of the cores. 
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I. INTRODUCTION 

Nanoalloys are bi- or multi-component metallic particles in the size range between 1-100 nrr>^. 
These nanoparticles are of great interest for both basic science and technological applications 
to magnetism, catalysis and optics. This wealth of possible applications stem from the very high 
degree of tunability of the physical and chemical properties of these systems, which is a direct con- 
sequence of the great variety of morphologies that nanoalloys can assume. Nanoalloys can present 
a variety of geometric shapes, ranging from crystalline structures, i.e. fragments of bulk crystals, 
to non-crystalline structures, the most common being icosahedra, decahedra and polyicosahedra^. 
But variety is not limited to the geometric shape, since the atomic arrangement in nanoalloys can 
form qualitatively different chemical ordering patterns. Chemical ordering patterns can be inter- 
mixed (randomly or phase ordered), core-shell and multishell. Also phase separated patterns have 
been found, such as the so-called as Janus particles^'"*, and the ball-and-cup configurations^-. This 
variety of geometries and patterns reflects the complexity of nanoalloy energy landscapes^, which 
renders the theoretical prediction of most stable structures quite challenging. The determination 
of the most stable structures is the starting point for calculating the properties of nanoalloys, that 
in general depend both on geometry and on chemical ordering. The interplay of geometry and 
chemical ordering is far from trivial, and it is the subject of active investigation. 

Silver- and gold-containing nanoalloys are presently attracting notable interest. In particular, 
the cases of poorly miscible pairs, such as AgCu, AgCo, AgNi, and AuCo, have been investigated, 
both by experiments and simulations, by several groups that have explored a variety of properties 
and potential applications^i^-^— . AgCo, AgNi, and AuCo nanoalloys are studied for their inter- 
esting magneto-optical properties which stem from the fact that they are composed by a ferro- (Co 
or Ni) and a non-ferromagnetic metal (Ag or Au), the latter presenting a sharp surface plasmon 
resonance (SPR). The frequency of the SPR in these nanoalloys depends on composition and or- 
dering pattern, so that it can be tuned to a good extent ''^. Applications to catalysis have been also 
proposed. For example, AgCo nanoalloys have been employed in the oxygen reduction in alkaline 
media"*^, and AgCu clusters are proposed as more efficient catalysts than the corresponding pure 
clusters for oxygen reduction reactions^. AgCu nanoalloys find application in the fabrication of 
Pb-free solder interconnects, that can melt at low temperatures^. Core-shell Au-Co nanoparticles 
have been proposed also for biomedical applications, such as thermal ablation therapies and drug 
delivery, for which the magnetic Co core is coated by a shell of Au, which is a biocompatible 



material that can be easily functionalized^"^. 



All these systems present some common features. Besides being poorly miscible in bulk 
phases"^^, these systems have the element of larger atomic radius (Ag or Au) presenting smaller 
cohesive and surface energies than the other element. These features are all in favour of Ag or Au 
surface segregation in these nanoparticles, a behaviour that has been confirmed, both experimen- 
tally and computationally, by several studies^i^"^— i^-i^^-^^^— i^. However, the fact that the 
cluster surface is expected to contain mostly Ag or Au does not determine completely chemical 
ordering. As we will see in the following, the latter is strongly influenced by a complex interplay 
with the nanoparticle geometrical shape. 



Generally speaking, even though several interesting results have been produced, the relation be- 
tween geometry and chemical ordering in determining the most stable nanoparticle configuration 
is still to be understood in all these nanoalloys. In particular, it is still unclear how the geometric 
shapes influence the preferential chemical ordering. In this paper we propose a systematic com- 
putational investigation of this topic, for AgCu, AgNi, AgCo, and AuCo nanoalloys. We consider 
Ag- or Au-rich compositions, that, as we will see, present the most complex and interesting in- 
terplay between shapes and patterns. From the computational point of view, most optimization 
studies about these nanoalloys have treated so far small cluster sizes, below 100 atoms, whose ex- 
perimental observation at the electron microscope is quite challenging. There are very few global 
optimization studies for sizes above 100 atoms for these systems^--^. These studies are however 
focused on a specific structural motif (the chiral icosahedron) that is relevant for intermediate and 
Ag- or Au-poor compositions. Here we consider clusters of sizes up to more than 10^ atoms (di- 
ameters up to 3-4 nm), that can be more easily observed in electron microscopy experiments. This 
is especially interesting because techniques such as STEM with Z-contrast imaging^ and energy- 
filtered TEM^ allow to reveal the internal structures of binary nanoparticles. We show that shape 
and placement of the inner core is clearly correlated with the overall nanoparticle shape, cen- 
tered cores being possible in icosahedra but not in decahedra and fee nanocrystals. In icosahedra, 
the core is centered until it reaches a critical size at which a morphological instability develops. 
Finally, we show that multicenter cores are found in polyicosahedra. 



II. MODEL AND METHODS 

From now on, the label S will refer to shell atoms (i.e. Ag or Au) and the label C to core atoms 
(Cu, Ni or Co). 

The nanoparticles are described within an atomistic model developed within the second- 
moment approximation to the tight-binding model (SMATB potential), known also as Gupta 
or RGL potential^^. In this model, the potential energy of the system depends on the relative 
distances between atoms r^j and it is written as the sum of single-atom contributions Ei'. 
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The parameters pij, qij, Aij, ^ij, r^ij depend indeed on the atomic species of the pair only. There- 
fore pij can be either pss or pcc or psc, and the same holds for all other parameters. As for 
i^oij, ^occ and tqss are the equilibrium distances in the bulk pure crystals at zero temperature, 
and tqsc = {rocc + ^055) /2. r^ is an appropriate cutoff distance so that the sum has contribu- 
tions from all neighbors of atom i that are within r^. In this work we choose to put r^ equal to 
the second-neighbor distance in the respective bulk solids for SS and CC pairs, whereas for SC 
interactions the cutoff distance is the arithmetic average. The potential is then linked to zero at 
the third-neighbor distance by a polynomial function in such a way that the resulting function is 
continuous with continuous derivatives. The parameters of the potential can be found in^^ for 
AgCu and AgNi, in^ for AgCo and in^ for AuCo. The interaction potential has been favorably 
checked against density-functional calculations^^^. 

We note that all systems present a considerable lattice mismatch between the metals. If we 
measure the mismatch by the quantity {r^ss — rocc)/ross we obtain 0.118 for AgCu, 0.138 for 
AgNi, 0.135 for AgCo and 0.131 for AuCo. 

In the following we will calculate the local pressure. For atom i, the local pressure Pi is given 

by 

dEi 

^' = -# <^' 

where Ei is the atomic energy and Vi is the atomic volume. Pi is expressed as a function of the 
atomic coordinates as 

(3) 
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where Vj is the atomic volume in the bulk solid and 



D{z) = (4) 



The functions (pij and V'ij are given by 
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if Tij is within the second-neighbor distance, by a polynomial link between the second- and third- 
neighbor distance and zero elsewhere. 

The lowest energy structures are searched for by a basin-hopping procedure^ with exchange 
moves only. In an exchange move, the positions of two atoms of different species are swapped. 
The move is followed by local relaxation to reach the closest local minimum. Both random and 
tailored exchanges are used^^. In tailored exchanges, atoms are chosen with different probabilities 
according to their local environment. Usually tailored exchanges allow finding chemical ordering 
arrangements with lower energy than random exchanges. Random exchanges become not efficient 
at all when dealing with the nanoparticle sizes considered in this work. It has been checked 
however that at smaller size, both random and tailored exchanges lead to the same results, the 
latter being faster in reaching the lowest-energy structures. For each size and composition, at least 
five independent simulations of 10^ steps each have been made. 



Ill, RESULTS AND DISCUSSION 

The main geometrical motifs that have been singled out in AgCu, AgCo, AgNi and AuCo clus- 
ters are fee nanocrystals, icosahedra, decahedra, and polyicosahedr a'-^'^^i^^-^^i^° . Representative 
clusters of these geometrical motifs are shown in Fig. [1] In the following, we consider each geo- 
metrical motif and we look for the most energetically stable chemical ordering within that motif, 
for a series of selected compositions. Our global optimization of the chemical ordering is made 
by means of the algorithms described in the Methods section, within the atomistic potential model 
whose form and parameters are given there. 
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FIG. 1. (A) Fee tmneated oetahedron of 1289 atoms. (B) Maekay ieosahedron of 1415 atoms. (C) Marks 
deeahedron of 1228 atoms. (D) Polyieosahedron of 419 atoms. It is made of five interpenetrating ieosahedra 
of 147 atoms. 

A. Fee struetures 



This motif has been studied in a series of computational works^i^i^ showing that in AgCu, AgNi 
and AgCo the cores assume preferentially off-center positions in clusters in the size range of 2 nm. 
This has been recently confirmed experimentally for AgCu^ in the case of larger nanoparticles, 
with diameters in the range 5-10 nm. As we show in the following, our calculations reproduce 
this behavior for sizes close to the experimental ones, and for all systems including AuCo. We 
consider specifically a truncated octahedron of 1289 (diameter ~4 nm) atoms and change the 
proportion of C atoms, explicitly considering the cases of 1 , 2 , 5, 20 and 200 C atoms. 

For a single C atoms, the most favorable positions are susburface sites in all cases (see the data 
reported in Table H referring to the impurity sites of Fig. |2]). This confirms what was already 
found for AgCu and AgNi fee clusters by means of atomistic calculations^ and for AgCu clusters 
by DFT calculations^. The most favorable position is always the subsurface site below a cluster 
vertex, with energy gains with respect to central inner sites that vary from a minimum of 0.057 



eV for AuCo to a maximum of 0.170 eV for AgCo. Surface sites are disfavored compared to 
subsurface and inner sites, especially in AuCo. The energetic cost of bringing an impurity atom 
from an inner site to the center of a (100) facet varies from 0.05 eV in AgCo, to 0.09 eV in AgCu, 
to about 0.15 eV in AgNi and to more than 0.3 eV in AuCo. Therefore in AgCo the energetic cost 
for placing atoms in the cluster surface is especially small. This has some important consequences 
that we will analyze in the following. 
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FIG. 2. Single-impurity sites in a truncated octahedron of 1289 atoms considered in Table U The vertices 
of the truncated octahedron are represented by midsize spheres, whereas other atoms than the impurity are 
shown as small spheres. 

For inner and subsurface sites, in which the nearest-neighbor bond numbers are always the 
same, the results in Table |I] show that the energetics is in general well correlated with the local 
pressure P acting on the C atom. P is always negative, indicating the expected tensile strain for 
impurities that are hosted in a matrix of bigger atoms. The most favorable site has the lowest 
absolute value of P in all systems, indicating that the lowering of strain is an important driving 
force in selecting the most favorable sites. For what concerns other inner sites, the correlation 
between energy and low absolute pressure is perfect for AgCu, AgNi and AgCo, while there is 
some discrepancy in AuCo for what concerns the pressures of the second- and third-best sites. 

For two C atoms, the most favourable configuration is always a dimer with one atom placed at 
a subsurface vertex position. Configurations with the two C atoms in dissociated configurations 
(for example, the C atoms placed at two different subsurface vertices) are always slightly higher 
in energy, by a minimum of 0.01 1 eV in AuCo to a maximum of 0.037 eV in AgNi. 

The case of five C atoms is quite interesting because of the possible appearance of a C atom at 
the cluster surface. In fact, for AgCu, AgNi and AgCo the lowest-energy configuration contains 
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surface, center of a (100) face 


0.087 


-6.59 


0.146 


-6.82 


0.050 


-5.27 


0.324 


-12.41 


surface, center of a (1 1 1) face 


0.083 


-7.48 


0.141 


-9.17 


0.050 


-7.51 


0.285 


-12.33 


inside the cluster, central site 


0.000 


-10.05 


0.000 


-13.76 


0.000 


-11.22 


0.000 


-15.40 


1 subsurface, below a vertex site 


-0.068 


-6.46 


-0.138 


-10.21 


-0.170 


-6.88 


-0.057 


-9.51 


2 subsurface, below the center of a (100) face 


-0.029 


-8.29 


-0.062 


-11.98 


-0.081 


-9.18 


-0.001 


-12.81 


3 subsurface, below the edge of a (100) face 


-0.049 


-7.38 


-0.101 


-11.10 


-0.126 


-7.97 


-0.030 


-11.19 


4 subsurface, below the edge of a (1 1 1) face 


-0.042 


-8.12 


-0.089 


-11.85 


-0.100 


-9.04 


-0.050 


-12.54 


5 subsurface, below the center of a (11 1) face 


-0.026 


-9.01 


-0.051 


-12.78 


-0.056 


-10.23 


-0.033 


-13.78 



TABLE I. Energy E and local pressure P of a single impurity in some different positions in a truncated 
octahedron of 1289 atoms, after local relaxation. The energy in the central position inside the cluster is 
set to zero in all systems. Subsurface positions are numbered as in Fig. |2l Energies ai^e given in eV, local 
pressures in GPa. 



a square-basis pyramid of C atoms which is placed as shown in Fig. [3tA)-(B). The vertex of 
the pyramid appears at the cluster surface, on a (100) facet. The origin of this effect can be 
qualitatively explained as follows. Let us compare this pyramid with the inverted pyramid in 
which the vertex atom is exchanged with an Ag atom and placed inside the cluster, as in Fig. 
[IC)-(D). 

Therefore, the transformation from the inverted pyramid to the pyramid with the vertex at the 
surface amounts in total to bringing a C atom (the vertex of the pyramid) from an inner site to 
the (100) facet, which should have an energy cost in the range 0.10-0.15 eV for AgCu, AgNi and 
AgCo (see the discussion about a single C atom). However, if one calculates the local pressure P 
on each atom of the pyramids, one obtains that its absolute values are well above 5 GPa for the 
inverted pyramid, while they drop down of about a factor 5 in the pyramid with the vertex at the 
surface. This better local relaxation overcompensates for the energy cost associated to bringing 
the C atom at the surface. The situation is different in AuCo, where this energy cost is much 
higher (about 0.3 eV) so that a better relaxation is not sufficient. In AuCo, the lowest-energy 
configuration has no C atoms at the cluster surface. 

For 20 C atoms the core is compact and off-center, with most of its atoms in the subsurface 
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FIG. 3. (A) and (B): side and top view of the lowest-energy structure of a cluster of 5 C atoms in a truncated 
octahedron of 1289 atoms, in the cases of AgCu, AgNi and AgCo. One C atom of the 5-atom pyramid 
appeai^s at the cluster surface, in a (100) facet. (C) and (D): side and top view of the 5-atom pyramid 
obtained by putting its vertex inside the cluster, so that no C atom appears at the surface. In (A) and (C) S 
atoms are shown as smaller spheres 



layer (see Fig. S]). In AgCu and AgNi there are two atoms in the surface layer, whereas there are 
8 in AgCo, and none in AuCo. The results found for AgCo is related to the low energetic cost 
related to placing Co atoms at the cluster surface. 

In Fig. |4]we show also the optimized structures of truncated octahedra of 1289 atoms containing 
cores of 200 atoms. Again, the cores are clearly off-center. In all cases, part of the core is covered 
by a single Ag or Au monolayer. Even though all systems present qualitatively similar behaviors, 
some differences can be noted. AgCo and AgNi present more compact cores, followed by AuCo 
and AgCu, which is the least compact. Moreover, in AgNi, and even more in AgCo, some C atoms 
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FIG. 4. Optimized structures of fee truneated oetahedra of 1289 atoms, with composition S1269C20 (top 
row) and S1089C200 (bottom row). S atoms are shown as small spheres so that C atoms (bigger spheres) are 
visible. 








FIG. 5. Left: unreconstructed (100) facet of the 1289-atom AuCo truncated octahedron. Right: recon- 
structed (100) facet with Au atoms covering the underlying Co core. 

appear at the cluster surface, especially in (001) facets, confirming the trends shown for smaller 
cores. These findings are in qualitative agreement with the fact that AgCo and AgNi have the 
strongest tendency against mixing^, and a large difference in cohesive energy with Ag. Therefore 
the Co and Ni cores in Ag try to maximize the number of C-C nearest-neighbor bonds. We note 
however that the number of Co or Ni atoms at the surface remains comparatively small (~10-15), 
because their appearance is mostly on the small (100) facets. This issue however deserves further 
investigation by fully unconstrained global optimization in which both cluster geometry and shape 
are optimized together^", so that shell atoms are allowed to better find their optimal configuration. 
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Work is in progress on this point. 

An interesting effect appears on the (100) facet of AuCo nanoparticle which is on the side 
where the core is placed. The Au atoms cover the Co core with a single layer which undergoes a 
reconstruction (see Fig. ID of the same type which occurs on Au(lOO) bulk crystal surfaces^^. 
This reconstruction leads to the formation of a more compact surface layer, with a hexagonal 
pattem instead of a square pattern. At variance with the case of infinite surfaces, we find this re- 
construction to occur only on the facets that cover the subsurface Co core, but not on other facets. 
The size mismatch between Au and Co atoms favors indeed Au atoms adopting more compact 
configurations, as those obtained by means of this reconstruction. This has been verified by arti- 
ficially increasing the size of Co atoms rocc while leaving all other parameters of the interaction 
potential unchanged. With no increase of tqcc, the configuration with unreconstructed facet is a 
local minimum, even though it is somewhat distorted and 1.8 eV higher in energy than the recon- 
structed facet. If tqcc is increased by 2%, the unreconstructed facet is still higher in energy (by 
0.6 eV), but it is not distorted. A reconstructed facet in which a single Au row is displaced be- 
comes slightly lower in energy than the configuration reported in the right panel of Fig. |5l For an 
increase of 5%, the unreconstructed facet becomes lower in energy than the reconstructed one, and 
finally, for an increase of 10% (at which C atoms are still smaller than S atoms), the reconstructed 
configuration is not even a locally stable minimum. The same tendency towards reconstruction is 
seen to a lesser extent in the AgCu case, confirming that an important driving force is indeed size 
mismatch, because Ag(lOO) bulk surfaces do not reconstruct. 

We note that in the reconstructed facet there are some Au atoms at the extremities of the rows 
that have rather low coordination, having five first neighbors only. This implies that displacing 
these atoms to other terrace positions has a small energy cost, of a few hundredths eV. It is thus 
expected that, with increasing temperature, these atoms may leave their original position to diffuse 
on the cluster surface. However, the facet should keep its reconstructed shape (with some shorter 
Au rows) because the unreconstructed geometry is much higher in energy. 

In all systems we have found a clear preference in favor of compact cores. However, Monte 
Carlo simulations of an atomistic model of AuPt were revealing the possibility of finding three- 
shell onionlike structures as the equilibrium ones^^. In our case, these structures, having non- 
compact aggregates of C atoms, are metastable^. In order to investigate the origin of the difference 
in behavior with respect to AuPt, we have verified whether it is possible to obtain onionlike struc- 
tures in AgCu by artificially changing the interaction parameters. For example, we have increased 
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FIG. 6. Lowest energy stmcture of AgnggCuioo for the case in which the mixed Ag-Cu interactions are 
artificially reinforced by increasing the parameter S^sc by 5%. Cu atoms fill subsurface positions below 
(100) facets instead of forming compact aggregates as those found the case in which the mixed interactions 
are not artificially reinforced. 



the parameter ^sc^ which rules the strength of the attractive part of the mixed interactions, by 
5%, finding that non-compact aggregates of Cu atoms fill the subsurface shell, especially below 
vertices, edges and (100) facets. These aggregates (see Fig. |6l reporting results for 100 C atoms) 
are very similar to those found for AuPt in Ref.^. However, putting a larger number of Cu atoms 
leads to filling both subsurface sites below (111) facets and inner sites. Therefore, the perfect 
three-shell structure, with the surface shell being completely of Ag atoms, the subsurface shell of 
Cu atoms and the inner part of Ag atoms, is not obtained as the lowest-energy configuration, but 
as a metastable one^. We note however that the surface shell is completely of Ag atoms. 
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B. Icosahedra 

Icosahedral core-shell nanoparticles present an even richer behavior, with the onset of a new 
kind of morphological instability. Before showing how this instability occurs, a few words about 
icosahedral structures in general are necessary. The Mackay icosahedron^"^ has a shell structure. 
An icosahedron with n shells (n > 1) has a number of atoms N = {lOn^ — 15n^ + lira — 3)/3, 
which gives the series of magic numbers 1, 13,55, 147,309,561,923, 1415 .... In the following we 
consider compositions for which there is a number of C atoms that corresponds to the completion 
of k shells, with k < n. 

At variance with fee and decahedra, the icosahedron is a structure whose central site is highly 
compressed. For monometallic clusters, the local pressure P on the central site is positive and 
quite large (several ten GPa) so that eliminating the central atom and leaving a vacancy at its 
place can be energetically favorable^^, if the icosahedron is sufficiently large. For the same reason, 
introducing a smaller atom at the center of the icosahedron instead of the vacancy, as in the case of 
Cu, Ni, Co in Ag or Au, allows a notable strain release^. The local pressure remains still positive 
(compressive strain), even if the impurity atom is hosted in a matrix of larger atoms, but strongly 
decreases compared to the case in which the center is occupied by an atom of the same species as 
the matrix. For the small impurity, the central site thus becomes by far the most favorable placed. 
This specific feature of the icosahedron will have important consequences for the structures of the 
cores, as we show below. 

As a first case we consider an icosahedron with fixed n = 8, i.e. with 1415 atoms in total, and 
increase the number of C atoms (we consider 1, 13, 55, 147, 309 C atoms, corresponding to k 
from 1 to 5, as shown in Fig. |7]). For k = I (not shown in the figure), the single C atom is always 
preferentially placed at the center of the structure as expected. For k = 2, the 13 C atoms are 
again placed at the nanoparticle center in all cases. However, when k increases above 2, a different 
behavior appears in AgNi. The 55 Ni atoms are not all contained in the inner shells with A; < 3, 
but some of them are now placed in the 4*^ shell, with some Ag atoms being correspondingly in 
the 3'"'^ shell. For AgCu, AgCo, and AuCo, the C atoms still form a perfectly centered icosahedral 
core. 

When k increases to 4 and 5, the morphological instability of the cores becomes more and 
more evident, with the development of irregular shapes that extend asymmetrically from the cen- 
ter towards the surface of the nanoparticle, filling some parts of the 6-th and 7-th shells. The 
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FIG. 7. Lowest-energy configurations of icosahedral nanoparticles of fixed size (number of shells n = 8, 
total number of atoms A^ = 1415) for increasing size of the core. From top to bottom: core of 13 atoms 
(number of shells k = 2), of 55 atoms {k = 3), of 147 atoms {k = 4), of 309 atoms (k = 5). S atoms are 
shown as small spheres so that C atoms (bigger spheres) are visible. 



asymmetric core however does not extend to the 8-th shell because C atoms gain in occupying 
subsurface positions, as we have seen for fee clusters^, so that none or very few C atoms appear 
at the cluster surface. The displacement of the core from the symmetric position is quantified in 
Fig. [8]by measuring the distance dc of the core geometric center from the central site of the icosa- 
hedron. The instability begins when dc shows a non-zero value. One can distinguish a weak form 
of morphological instability, when only sites of the (k + 1)*'^ shell are occupied by C atoms and 
dc has a small value, from a fully developed form, which is associated to much larger values of 
dc. 

We remark that centered symmetric core shapes still remain locally stable minima at least up to 
/c = 5 in all cases. However, with increasing k, the energy difference between the most favorable 
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FIG. 8. Displacement dc (in A) of the geometric center of the core with respect to the center of the 
icosahedron, as a function of the number of shells in the core k. Different symbols refer to different systems 
(+ to AgCu, X to AgNi, * to AgCo, and □ to AuCo). 





AgCu 


AgNi 


AgCo AuCo 


nk 


A^ 


A^ 


AE 


/\E 


8 1 


0.00 


0.00 


0.00 


0.00 


8 2 


0.00 


0.00 


0.00 


0.00 


8 3 


0.00 


-0.07 


0.00 


0.00 


8 4 


-0.52 


-2.06 


-0.89 


0.00 


8 5 


-4.34 


-14.85 


-13.55 


-7.05 


64 


-1.53 


- 


- 


-2.69 


7 4 


-0.85 


- 


- 


-0.40 



TABLE II. Energy differences /S.E (in eV) between the lowest energy structure and the structure with 
centered symmetric core. Negative /S.E indicate that asymmetric cores are more favorable. 



asymmetric core and the metastable centered symmetric core increases, becoming very large for 
k = h especially for AgNi and AgCo (see Table II]). 

These results show that there is a critical value of k, referred to as kc, at which the morpholog- 
ical instability starts, kc is system-dependent. For n = 8, /Cc = 3 for AgNi, kc = ^ for AgCu and 
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AgCo, kc = 5 for AuCo. 

Now we consider the possible dependence of kc on n. In Fig. |9]we fix /c = 4 and increase 
n from 6 to 8. The cores that are asymmetric for n = 6 become more and more centered and 
symmetric for n = 7 and 8, showing that increasing the shell thickness tends to stabilize the 
morphology of the cores. This shows for example that, in AuCo, /cc = 4 for n = 6 but it increases 
to kc = 5 for n = 8. 

What is the physical origin of this morphological instability? The instability has the same 
qualitative features in all systems, with some quantitative differences, developing earlier in AgNi, 
then in AgCu and AgCo and finally in AuCo, the latter system being able to support the largest 
centered symmetric cores. The onset of the instability in all these systems suggests that there 
might be an explanation of generic character, originating from some simple common feature, and 
being possibly applicable to an even wider class of systems. 

The key property that has to be considered turns out to be the lattice mismatch between the 
elements in the nanoparticle. This can be understood by quite simple considerations. In the bulk 
crystals, the nearest neighbor distances between S and C atoms are tqs and tqc, respectively. Now 
we change artificially roc in our model letting all other parameters of the model unchanged. This 
amounts to changing the size of C atoms only. Then we look for the value r^^ which minimizes 
the energy of the nanoparticle, and we compare it to the original rocr^. We repeat this procedure 
for icosahedra of 1415 atoms (n = 8) with centered symmetric icosahedral cores with different k, 
from 1 to 5. 

Let us consider for example the case of AgNi, for which r^s = 2.89 A and tqc = 2.49 A. 
When k = I the single Ni atom in the central site has an optimal size Tqq which is smaller than 
roc by 8.3%. This means that the strain release which is achieved by placing a Ni atom at the 
center of an Ag icosahedron would be even more effective if the size of the Ni atom would be 
smaller than its actual size. The situation begins to change for A; = 2. The 13 Ni atoms form a 
small central icosahedron. However, in this case r^^ is slightly larger than tqc, by 2.7%, so that 
the best strain release would be achieved if Ni atoms were slightly bigger. This tendency becomes 
more evident for A; = 3, for which Tqc is larger than roc by 5.7%. This means that the actual size 
of Ni atoms is indeed too small to support a centered icosahedral core of 55 atoms . This core is 
therefore too small, so that it may benefit by the substitution of some of its Ni atoms by larger Ag 
atoms, as shown in Fig. |71 second row. The behavior of AgCu and AgCo is of the same type, but 
the instability occurs at larger core sizes due to a smaller lattice mismatch between the constituent 
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elements. The same line of reasoning apply to AuCo too, even though the comparison with the 
other systems is less direct due to the fact S atoms are Au. 

The physical effect triggering the instability can be analyzed also by calculating the local pres- 
sure Pi acting on each atom, as given by Eqs. ^^. 

In Fig. [To] we report P^ for the inner 55 atoms for the 1415-atom icosahedron in the case of 
AgNi with centered symmetric icosahedral cores. For pure Ag (k = 0) Pi is strongly positive 
at the central site, which is therefore strongly compressed, and, on average decreases with the 
distance from the center—, being however positive for all atoms belonging to the first three shells. 
Substituting a single Ni impurity decreases Pi at the central site^, but the local pressure in other 
sites remains nearly unchanged. For k = 2 we have two shells of Ni atoms, and we begin to 
see some qualitative change in the behavior of Pi. In fact, the pressure of the 12 Ni atoms in the 
second shell becomes negative, indicating again that these atoms show some tendency to expand. 
Finally, for A; = 3, Pj becomes more strongly negative for the atoms of the third shell, especially 
for those that are at vertex sites of the 55-atom icosahedron. In this case it is favorable to substitute 
these vertex atoms with bigger atoms, thus obtaining a small energy gain. In the lowest-energy 
structures we have thus a few Ag atoms in the third shell in vertex position and, correspondingly, 
the same number of Ni atoms is now forming a small island in the fourth shell. 

This morphological instability is thus originating from the accumulation of strain as the number 
of shells in the core increases. The key point is that strain changes from a compression for small 
cores to a stretching for large cores. Since S atoms are of bigger size, this kind of strain can 
be somewhat released by substitution of C atoms in the inner shells. This leads to an energy 
gain, which is small for the first substitutions but becomes quite large as the number of C atoms 
increases, as shown in Table HIl 

The morphological instability of the cores therefore bears some common points with the 
Stranski-Krastanov (SK) instability occurring in heteroepitaxial thin film growth^. In heteroepi- 
taxial growth, atoms of species B are deposited on a substrate of species A. Between the two 
species a lattice mismatch is present. As the number of deposited layers increases, strain accu- 
mulates so that the film breaks up in its top part, the growth mode changing from layer-by-layer 
to three-dimensional, with the formation of mounds. The SK instability occurs when there is a 
(free) energy gain in starting the (n + 1)*'* layer before completing the n^'^ layer—. Therefore both 
SK instability and core instability in icosahedra share a common origin. There are however some 
differences. 
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FIG. 9. Lowest-energy configurations of icosahedral nanoparticles of fixed core size (number of shells 
k = 4, corresponding to 147 C atoms) for increasing number of S atoms (and size of the nanoparticle). 
From top to bottom: N = 561 (number of shells n = 6), N = 923 (n = 7), N = 1415 (n = 8). S atoms 
are shown as small spheres so that C atoms (bigger spheres) are visible. 
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FIG. 10. Local pressure Pi in an icosahedron of 1415 atoms, for its inner 55 atoms that are arranged with 
increasing distance from the center. Different symbols refer to different compositions, i.e. different number 
of Ni atoms. +: pure Ag cluster (k = 0). x: single Ni impurity at the central site {k = 1). *: centered 
icosahedron of 13 Ni atoms (k = 2). D: centered icosahedron of 55 Ni atoms (k = 3). The bottom pictures 
show the difference between the centered symmetric (left) and the optimal (right) cores in the case of 55 
Ni atoms. In the optimal core, some Ag atoms substitute the Ni atoms at vertex positions in the third shell. 
These Ni atoms are now placed in the fourth shell, which is shown by smaller spheres. 

In the SK instability, the strain corresponds to either compression or stretching, depending 
on the sign of the lattice mismatch. In the core instability, the strain changes sign, from the 
compression, which is specifically inherent to the icosahedral geometry, to the stretching which 
originates by the increases of the number of shells in the core. Therefore the core instability 
requires a specific sign of the lattice mismatch, with C atoms being smaller than S atoms. 

Moreover, in the SK instability, the strain is caused by the lattice mismatch at the bottom layer 
of the thin film, and it is released by breaking the thin film in its outer layers. In the core instability, 
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strain is originated by lattice mismatch at the outer shell of the core, and it is released by breaking 
these outer shells. We note that concerning the critical core size, we have determined a critical 
number of shells kc for its onset. More precisely, the instability starts from a number of C atoms 
Nc which is between the numbers of atoms necessary to complete the shell with kc — I and the 
shell with k^. 

We remark that the term morphological instability does not mean that centered cores are un- 
stable configurations. They are locally stable minima, whose energy is however much higher than 
the energy of the optimal configuration. The term morphological instability means that the cen- 
tered morphology does not persist as the lowest-energy morphology as the proportion of C atoms 
increases in the cluster, i.e. the term instability is referred to the overall core morphology and not 
to the local minimum of the centered core. 

We briefly discuss now what happens going towards C-rich compositions. For example, we 
have considered the case k = n — 1 in AgCu. For Cu-rich cases, it was previously shown that the 
Mackay structure of the shell is not at all the optimal one, anti-Mackay or chiral shells being much 
lower in energy^i^^. Essentially, the Mackay shell contains too many atoms to well adapt itself to a 
"small" substrate. Anti-Mackay or chiral shells contain fewer atoms so that they can adapt much 
better. If the number of Ag atoms is that of a Mackay shell, the optimal configuration is indeed 
the anti-Mackay or chiral shell with the Ag atoms in excess forming an island above it. Therefore 
it is now the shell which loses its symmetry, while the core can stay symmetric. However it is 
possible to build a Mackay shell for k = n — 1, and then optimize chemical ordering by exchanges 
only, thus keeping this geometric structure. This leads to an intermixing of Ag and Cu within the 
surface shell ( which is heavily strained if made of Ag only), with some Ag atoms correspondingly 
going inside the cluster. 

Finally consider icosahedra with a central atomic vacancy—, that might be more stable even in 
the bimetallic case especially at large sizes. We consider in details AgCu with n = 8, correspond- 
ing to a cluster of 1414 atoms. It turns out that the behavior is quite the same as in the case without 
the vacancy. For k = 2 and 3 (12 and 54 Cu atoms respectively) we find centered symmetric 
icosahedral core morphologies with a central atomic vacancy. From k = A on, the morphological 
instability takes place, as in the case of no vacancy. 
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C. Decahedra 

The most favorable type of decahedron for pure Ag and Au clusters of sizes in the range ~ 10^ 
atoms is expected to be the Marks decahedron^. Decahedral clusters have a single fivefold axis, 
which is shared by the five tetrahedra composing the structure. The axis terminates into two 
fivefold vertices. Most part of the decahedron surface is made of (1 1 l)-like facets, but there are 
five (lOO)-like open facets too. In the following we focus on a Marks truncated decahedron of 
1228 atoms. 
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FIG. 11. Single-impurity sites in a decahedron of 1228 atoms considered in Table Hill 





AgCu 


AgCu 


AgNi 


AgNi 


AgCo 


AgCo 


AuCo 


AuCo 


position 


E 


P 


E 


P 


E 


P 


E 


P 


1 subsurface, below a fivefold vertex 


-0.063 


-4.47 


-0.137 


-8.23 


-0.169 


-4.61 


-0.161 


-7.46 


2 subsurface, below the edge of a (100) facet 


0.013 


-6.49 


-0.025 


-10.16 


-0.068 


-7.13 


-0.041 


-9.54 


3 subsurface, below the edge of a (1 1 1) facet 


0.032 


-7.81 


0.012 


-11.44 


-0.012 


-8.46 


0.002 


-11.77 


4 subsurface, below the center of a (100) facet 


0.054 


-8,30 


0.052 


-12.00 


0.022 


-9.17 


0.099 


-12.78 


5 subsurface, below the center of a (1 1 1) facet 


0.060 


-9.22 


0.067 


-13.10 


0.051 


-10.49 


0.096 


-13.89 


6 central, on the fivefold axis 


0.000 


-7.94 


0.000 


-11.60 


0.000 


-8.78 


0.000 


-13.42 


7 central, nearest neighbor of the fivefold axis 


0.038 


-9.36 


0.055 


-13.13 


0.053 


-10.66 


0.063 


-14.87 



TABLE III. Energy E and local pressure P of a single impurity in some different sites in a decahedron of 
1228 atoms, after local relaxation. The sites are represented in Fig. [TT] The energy in the central position 
inside the cluster is set to zero in all systems. Energies are given in eV, local pressures in GPa. 

By analyzing the local pressure in a pure metal Ag or Au decahedron, it turns out that the sites 
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with the highest local positive pressure are subsurface sites below the fivefold vertices, which thus 
belong to the fivefold axis. Their local pressure is not so high as in the central site of the pure-metal 
icosahedron, but it is still notable (8 and 16 GPa in the Ag and Au decahedron, respectively). It 
is therefore not surprising that these are the most favorable sites for small impurities, as shown 
in Table Hill and Fig. [TTJ Other favorable sites are either along the fivefold axis or in subsurface 
positions, especially below (lOO)-like facets. 



Agi223CU5 Agi2i7CUji AgisisCUjs AUi20sCo 
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FIG. 12. Optimized structures of Marks decahedra of 1228 atoms and different compositions. Each cluster 
is shown in top and side view (first and second row, respectively). S atoms ai^e shown as small spheres so 
that C atoms (bigger spheres) ai^e visible. 

Favorable sites along the fivefold axis are arranged along a straight line, therefore it is not 
possible to build up compact aggregates by filling them only. On the other hand, this is possible 
by filling subsurface sites below (100) facets. This generates an interesting competition at small 
core sizes between the tendency to fill the most favorable sites and the tendency to form compact 
aggregates. This competition is especially evident in AgCu. For Agi226Cu2, the two subsurface 
sites below vertices are occupied. Increasing the number of Cu atoms to obtain Agi223Cu5, two 
small linear aggregates (a dimer and a trimer) are formed along the fivefold axis as shown in Fig. 
[T2I The fivefold axis is filled by Cu atoms up to composition Agi2i7Cuii, at which all inner sites of 
the axis are filled. Up to this composition, the tendency to fill the most favorable sites prevails over 
the tendency to form compact aggregates. However, a further increase of Cu causes the tendency 
in favor of compact aggregates to prevail. At composition Agi2i5Cui3, a transition to completely 
different core shapes, that are more compact and below (100) facets, finally takes place. In AgNi 
and AgCo, the behavior is of the same kind, but the transition to compact shapes occurs at even 
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FIG. 13. Optimized structures of Marks decahedra of 1228 atoms, with compositions SnygCso (top row) 
and S1028C200 (bottom row). S atoms are shown as small spheres so that C atoms (bigger spheres) are 
visible. 



smaller core sizes. 

On the other hand, AuCo shows a different behavior. In fact, the small compact cores do not 
form below a (lOO)-like facet, but in subsurface position below a fivefold vertex, as shown in Fig. 
[l2]forAui208Co2o. 

However, in all systems, large cores form below (lOO)-like facets as shown in Fig. [13] for 
cores of 50 and 200 atoms. These cores are clearly off-center and somewhat resemble those found 
in fee nanoparticles. The cores are slightly more compact for AgNi and AgCo. For AgCu and 
AuCo the most favorable core shapes of size 200 are quasi-symmetric with respect to a twin plane, 
resembling two incomplete tetrahedral units that share a few atoms on the twin plane. 



D. Polyicosahedra 

Polyicosahedra are made of several icosahedra that share some atoms. Core-shell polyicosahe- 
dra have been theoretically predicted^ as being very stable structures for AgNi and AgCu clusters 
for small sizes, below 50 atoms. For these sizes, the polyicosahedra are made of interpenetrating 
small icosahedra of 13 atoms, which in turn consist of a central atom (which is always either Cu or 
Ni) and of its first coordination shell, containing 12 atoms that can be in general of both species. 
However this is not the only way of building polyicosahedra. In fact, for larger sizes, one can 
consider larger icosahedra as elementary building blocks^. In the following, we will consider the 
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Mackay icosahedron^ of 147 atoms as the elementary building block, that allows to build up clus- 
ters of a few hundred atoms. In particular, we deal with a pentaicosahedron of 419 atoms (see Fig. 
[T])- If C atoms are few, this structure heavily distorts in the case of AuCo. Therefore we consider 
only AgCu, AgNi and AgCo in the following, unless otherwise specified. The results concerning 
the most stable structures are shown in Fig. [14] 





AgCu 


AgCu 


AgNi 


AgNi 


AgCo 


AgCo 


position 


E 


P 


E 


P 


E 


P 


1 icosahedral center 


0.000 


10.02 


0.000 


5.84 


0.000 


10.48 


2 segment joining icosahedral centers 


0.381 


-0.77 


0.494 


-4.86 


0.467 


-1.28 


3 segment joining icosahedral centers 


0.410 


-1.65 


0.527 


-5.69 


0.505 


-2.05 


4 neighbor of icosahedral center 


0.477 


-3.36 


0.621 


-8.07 


0.574 


-4.66 


5 center of the triangle of icosahedral centers 


0.526 


-4.74 


0.681 


-8.60 


0.648 


-5.36 


6 neighbor of icosahedral center 


0.569 


-6.20 


0.735 


-10.11 


0.704 


-6.97 


7 subsurface position 


0.583 


-3.41 


0.687 


-7.76 


0.562 


-4.20 



TABLE IV. Energy E and local pressure P of a single impurity in some different positions in a polyicosa- 
hedron of 419 atoms, atoms, after local relaxation. Energies are given in eV, local pressures in GPa. Site 
numbers refer to Fig. [15] 



The pentaicosahedron has five icosahedral central sites, that are separated from each other by 
two sites in between. Instead of forming a single aggregate, C atoms begin to populate the dif- 
ferent centers, forming very small aggregates, until there are atoms enough to merge in a single 
connected aggregate. As can be seen in the first row of Fig. [14] the way by which the different 
small aggregates merge is system-dependent. In AgCu, a non-compact network-like aggregate 
forms, with Cu chains joining the icosahedral centers. In AgNi and AgCo, a more compact ag- 
gregate of two layers forms, comprising three icosahedral centers, plus two isolated atoms in the 
remaining icosahedral centers. For 21 C atoms (second row) AgCu continues the network-like 
grown sequence, whereas AgNi and AgCo form a 20-atom tetrahedron (comprising four icosahe- 
dral centers) plus an isolated atom in the fifth center. Finally, for 30 C atoms, all systems form the 
same structure, a single connected aggregate with shape of a double tetrahedron whose five ver- 
tices are the icosahedral centers. Increasing the number of C atoms further, the growth sequence 
for AgCu and AgNi is more symmetric, with a chain of atoms developing around the common 
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FIG. 14. Optimized structures of polyicosahedra of 419 atoms. From top to bottom, the compositions are: 
S401C18, S398C21, S389C30, S365C54, and S323C96. S atoms are shown as small spheres so that C atoms 
(bigger spheres) are visible. 
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FIG. 15. Single-impurity sites in a polyicosahedron of 419 atoms considered in Table JV] Site 1 is an 
icosahedral center, sites 2 and 3 are on segments joining icosahedral centers, sites 4 and 6 are neighbors of 
an icosahedral center, site 5 is in the middle of the triangle whose vertices are icosahedral centers, site 7 is 
subsurface. The four icosahedral centers besides site 1 are indicated by mid-size spheres. Other sites are 
indicated by small spheres. 

basis of the two tetrahedra forming the double tetrahedron. For AgCo, the structures are less 
symmetric. However, in all cases, small 13-atom icosahedra of the C species form around the 
icosahedral centers. Finally, for 96 C atoms the growth sequence is completed by a cluster which 
preserves the symmetry of the original pentaicosahedron. Around each icosahedral center there is 
a small 13-atom icosahedron of C atoms. This last structure is stable also for AuCo. 

The growth sequences of the cores can be understood by considering site energies of single 
impurities and their local pressures, as reported in Table FVl which refers to the sites shown in Fig. 
[T5l The five icosahedral are by far the preferred sites. C atoms in these sites present a positive 
pressure, which is however by far smaller than the pressure that would act on an Ag atom in the 
same sites. Therefore, as in simple icosahedra, placing the impurity in a icosahedral center causes 
a notable strain release and the stabilization of the structure. The second-best sites for a single 
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impurity are then those on the segments joining the icosahedral centers. This is the driving force 
towards the formation of network-like aggregates. This driving force is counterbalanced by the 
tendency to form compact aggregates, which originates from the fact that the C-C bonds are the 
strongest ones. In AgCo and AgNi, where C-C bonds are by far stronger, the tendency towards 
more compact structures prevails. In AgCu, C-C bonds are less strong, therefore the driving force 
towards the network-like structures prevails up to about 20 Cu atoms. We note that low absolute 
values of the local pressures are well correlated with favorable site energies. 

IV. CONCLUSIONS 

Our results clearly show that the optimal shape and placement of the core in core-shell nanopar- 
ticles are strictly correlated to the overall geometry of the nanoparticle itself. However some com- 
mon features of general character can be singled out. The first point is that off-center asymmetric 
cores are the rule as lowest-energy structures, not the exception. In fact, we have seen that these 
cores are dominant in crystalline and decahedral motifs. Centered cores can be obtained in icosa- 
hedral nanoparticles, but only for core sizes that are not too large. In fact, our calculations show 
that there is a critical size for the onset of a morphological instability which drives the core to an 
asymmetric shape extending towards the nanoparticle surface. In polyicosahedral nanoparticles 
(and to a much lesser extent in decahedra), multi-center core structures are possible. Even in these 
cases the cores are often not symmetrically placed with respect to the center of the whole structure. 

In all cases in which off-center cores are prevailing, centered cores are still possible as 
metastable configurations, whose lifetime can drastically change depending on temperature. 

The physical factors causing the asymmetric placement of the cores and the onset of the mor- 
phological instability of the cores in icosahedra originate from very simple features of generic 
character, essentially from the strain originated from lattice mismatch. The morphological insta- 
bility of icosahedral cores bears some analogy with a very general phenomenon taking place in 
thin film heteroepitaxial growth, i.e. the Stranski-Krastanov instability. For this reason, we believe 
that the same effects should appear in a whole class of nanoparticles sharing the same features, i.e. 
poor miscibility, lattice mismatch and tendency of bigger atoms to segregate to the surface. For 
example, systems as AuNi, AuRh and AuPt are expected to present this kind of behavior. 

Even though our calculations are finding low-energy structures and therefore are related to the 
low-temperature behavior of the systems, we believe that the asymmetric cores should be present 
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also in a significant temperature range. In fact, energy differences between centered cores and 
off-center cores become readily very large, so that we expect that the effect of temperature may 
change the details of the core shape, however without changing its qualitative features. On the 
other hand, the melting behavior of icosahedra with small centered cores compared to the melting 
of icosahedra with larger off-center cores is an interesting subject which deserves investigation. 
Also the magnetic properties in Co- or Ni-containing nanoparticles with off-center cores can be of 
interest. 

Finally we note that our results are relevant to the design of experimental procedures to achieve 
effective coatings of the cores in systems in which lattice mismatch is present. Our results show 
indeed that achieving effective coatings may be a quite difficult task. In fact, in the lowest energy 
configuration, the core very often tends to approach the nanoparticle surface, so that part of the core 
itself is covered by a very thin layer of shell atoms. A few core atoms may even appear at the cluster 
surface in some cases. Centered cores can be experimentally grown, but these are metastable 
structures, which therefore may present aging problems due to this driving force towards more 
asymmetric shapes. 
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